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Abstract 

Using the volume coupling version of the cloudy bag model, the quark-meson coupling model 
is extended to study the role of pion held and the properties of nuclear matter. The extended 
model includes the effect of gluon exchange as well as the pion-cloud effect, and provides a good 
description of the nuclear matter properties. The relationship between the extended model and 
the EFT approach to nuclear matter is also discussed. 
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The quark-meson coupling (QMC) model [l| can be considered as an extension of Quan- 
tum Hadrodynamics (QHD) to include the effect of the internal structure of a nucleon in 
matter. The model describes a nuclear system by non-overlapping MIT bags, in which the 
confined quarks interact through the self-consistent exchange of isoscalar, scalar (a) and 
vector (uj) mesons. In the past few decades, it has been extensively developed and applied 
to various nuclear phenomena with tremendous success 

On the other hand, a major breakthrough occurred in the problem of nucleon-nucleon 
(NN) force by introducing the concept of an effective field theory (EFT) [2]. The QCD 
Lagrangian for massless up and down quarks is chirally symmetric, and the axial symmetry is 
spontaneously broken. This implies the existence of the massless Nambu-Goldstone bosons, 
namely the pions. The non-zero pion mass is then a consequence of the fact that the light 
quark has a small mass. Thus, one arrives at a low-energy scenario that pions and nucleons 
(and possibly deltas) interact via a force governed by spontaneously broken, approximate 
chiral symmetry. 

In EFT, the degrees of freedom of quarks and gluons (including heavy mesons and nu- 
cleon resonances) should be integrated out (or cutoff), because a probe with wavelength A 
is insensitive to details of structure at distances much smaller than it [3j. Instead, it is 
necessary to add local (contact) interactions with low-energy constants (LECs) to the La- 
grangian to mimic the effect of the true short distance physics. The LECs are determined 
empirically from fits to ttN and/or NN scattering data, and vary with the momentum cutoff 
(~ A -1 ) accounting for quantum fluctuations excluded by the cutoff [3[. Up to N 3 LO, EFT 
can provide the peripheral NN scattering data (below about 250 MeV lab. energy) very 
accurately [l]. 

Recently, EFT has been intensively applied to the problem of nuclear matter. In addition 
to the usual (small momentum) expansion in the free NN or 7rN scattering, physical observ- 
ables in matter are expanded in terms of the Fermi momentum fcp, which is also a relevant, 
small scale. Such density dependence arises from the Pauli blocking effect in matter, i.e., the 
medium insertion including the step-function, 6(kp— \p\), in the nucleon propagator. Then, 
the strengths of the LECs are fine-tuned so as to reproduce the nuclear matter properties jlj] . 

If the internal structure of the nucleon were completely frozen in a nuclear medium or 
the same as that in free space, it might be sufficient to consider the density dependence 
solely stemming from the Pauli blocking effect. However, if the in- medium nucleon were 
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metamorphosed depending on the nuclear density ps, the situation may be different. In fact, 
the evidence for the medium modification of nucleon structure was observed in polarization 
transfer measurement in the quasi-elastic (e, e'p) reaction at the Thomas Jefferson national 
accelerator facility, and the result supports the prediction of the QMC model [6[. It also 
seems vital to consider the internal structure change of the nucleon to understand the nuclear 
EMC effect pj. 

The QMC model can describe the medium modification of the nucleon structure through 
the quark model, and predict the density (or mean scalar-field) dependence of physical 
quantities Expanding such modification in terms of kp and comparing with the values of 
LECs given in the EFT approach, 1 it may be possible to study whether the internal structure 
change of a nucleon indeed shows up in matter, since the LECs involve all information on 
the short distance physics. 

To carry out such a complicated investigation, as a first step, we need to develop a new 
version of the QMC model for nuclear matter, where the structure of nucleon (and delta) is 
treated based on chiral symmetry. We attempt this in the present study using the volume 
coupling version of the cloudy bag model (CBM), which incorporates major results of the 



current algebra for low energy 7rN scattering 



9] 



The Lagrangian density for the volume coupling version of the CBM in flavor SU(2) is 



given by 
Ccbm 



9] 



6 v -~W5 s + ±(dJY + C xB , (1) 



with ip the quark field, $ the pion field, 4> = • 0) 1/2 , = }/\$\, U (= 93 MeV) the 
pion decay constant, B the bag constant, 6y the step function for the bag, 5s the surface 
5-function, Dj= (9^)#+/ w sin(0//,)^0 and = d^- |[cos(0/A) -l]r- (# x d»$)ij). 
The last term includes the quark mass, m, which explicitly breaks chiral symmetry, and the 
pion mass, m 7r (= 138 MeV): C xB = -rmpe'"'^ 5 ^ ^ipOy - \m* n $' 1 . 

As in the CBM, we linearize the pion field and keep 0(1/ f^) (the convergence properties 



1 In fact, from the poin of view of the quark and gluon degrees of freedom, the QMC model can explain 
the values of the coefficients appearing in the familiar (contact) Skyrme force in conventional nuclear 
physics 
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of the CBM were given in Ref. |9|). The Lagrangian density then reads 

CcBM '- 



— \ TYl -* 1 

i> | hixd" - m + i— 757* • + ^7m75^ • (<9 M 0) ^ V - B 



By - -"0#S 



(2) 



Here the pion field interacts with the quark through both the pseudovector (pv) and pseu- 
doscalar (ps) couplings. The strength of the ps coupling is 0(m/f n ), which explicitly shows 
the breaking scale of chiral symmetry. 

We introduce the gluon field as well. The resulting Lagrangian density is thus given by 



■>CBM 



BAG 



where 



£bag 

C-int = $ 



int ■> 



1 



i-f^ - m)tfj -B]&v- -^8 S , 



Tfl — » I — * Q —* —* 



(3) 

(4) 
(5) 



with A the SU(3) generators and g the quark-gluon coupling constant. The free pion field 
and the kinetic energy of the gluon field, A^ 1 , are, respectively, described by C n and C g . 

We firstly calculate the second-order energy correction to the nucleon or delta mass. 
The energy shift of a multi-quark, ground state, |0), due to the interaction is given by the 
Hubbard's prescription 



E-En 



(°iEht^/^)^i--- / 

m=l ' J J 



d^x m T\H int {xx i 



(6) 



where TLmt is the interaction Hamiltonian density. The energy shift is then given as E® = 
Edr + E n d, where the first term is the direct contribution and the second one is the non-direct 
contribution. (See Eqs. lfTTl) and (fl2l) later.) 

The noninteracting, quark green function is given by iG°(r,r') = {0\T[if)(r)il)(r')} |0), and 



it can be separated into two pieces: G°(r,r') = f 



du> iu>(t— t') \(~iQ 



The first term is the usual Feynman propagator in a spherical cavity (bag) and the second 



one describes the occupied, multi-quark ground state 



io| 



G° D (r,r',u) = U n (r)U n (r')2m6(u-E n ), 



(7) 



ii<«j 
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where U n is the positive energy state with a complete set of quark quantum numbers n{= 
{uKfifiific}) including isospin /ij and color fi c {rip specifies the quantum numbers at the 
Fermi surface in a hadron). 

Here we restrict the expansion of the quark propagator to the ground state, i.e., v — and 
k = — 1. Such a truncation may be considered as a regularization of the quark propagator, 
where in flavor SU(2) the intermediate baryon states in loop diagrams are restricted to the 



nucleon and delta This is consistent with the idea of the CBM. Thus, we let n label 
the spin, isospin and color {/i/ij/x c }. 

The pion propagator is defined by iA a b(r, r') = (O|T[0 a (r)0b(r')]|O) = i<5 a {,A(r, r'), where 
(a, b) specifies the isospin. It is then given by the multipole expansion 

A(r,r') =J2j ^e-^-^A e (ry,u;)YUf)Y; m (n- (8) 



The gluon propagator can be calculated in the Coulomb gauge 2 

iD^(r,r f ) = i8 cd J ^y~ tk T X) = iSoaDooM, (9) 

where (D^) represents the Coulomb (transverse) propagator, (z,j) = 1,2,3 and (c,d) 
specifies the color. 

Using these propagators, the direct (or Hartree) contribution between two quarks is 
calculated by 

E dr = ~\JJ S(t 1 )d' l r 1 d' l r 2 e(r 1 ,r 2 )tT[r 1 t 1 G D (ri,r 1 )] ■ tr[r 2 ^G° (r 2 ,r 2 )], (11) 

where Ti^ represents the vertex (ps, pv or gluon) including the coupling constant (the 
Lorentz index is suppressed here), 9 stands for the pion or gluon propagator and t = r/2 
or A/2. The non-direct contribution is given by 

E nd = l - Jf 5{t l )dyd i r 2 Q{r 1) r 2 )ti[Y l i l G\r 1 ,r 2 )-Y 2 t 2 G\r 2 ,r 1 )}. (12) 

This consists of the exchange contribution, E ex , which is evaluated with G° = G° D in Eq. fTl2l . 
and the self-energy contribution, E qs f, where the pion is emitted and absorbed by the same 



It can be shown that the result does not depend on the choice of the gauge 11 1. 
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quark. E qs f is calculated by replacing one of G° D 's with G° F in E ex . The case with simul- 
taneously G° = G° F is removed, because it is the vacuum energy. The Hartree-Fock (HF) 
contribution between two quarks is defined by E qHF = E dr + E ex . Note that these energy 
shifts correspond to the case where a closure approximation to the intermediate states is 
taken in the CBM, and that such an approximation is reasonable. 
For example, the HF result for the pv coupling, E^ F , is given by 



jat] HF N* r 2 r 2 r d tt4 

E aHF = - ln .,. Vf n^n / d PlPl / d P2p2 



12^(UR)mJ "™J "™J Q t 2 + y 



A(p 1 )A(p 2 )j (tp 1 )j (tp 2 ) 



+ - {A{p l )B{p 2 )j {t Pl ) (j (tp 2 ) - 2 ]2 {tp 2 )) + A(p 2 )B( Pl )j {tp 2 ) (j'o(tpi) - 2j 2 (t Pl ))} 



+ l -B{ Pl )B{p 2 ) ( 3o (t Pl ) - 2j 2 (t Pl )) ( 3o (tp 2 ) - 2 ]2 {tp 2 )) 



(13) 



with the normalization constant for the quark wave function, R the bag radius, x the 
lowest quark eigenvalue, y = m n R/x, A(p) = jg(p) — /3 2 jf(p), B(p) = 2f3 2 jf(p), ft 



xj [a + 5), a 2 = x 2 + 5 2 and 5 = rriR. Here the spin-isospin matrix element is given by [10] 



[^}hf= (A^')-^V^) = ^-S(S + 1)- 1(1 + 1), (14) 

where the index (i, i') runs over the spin and isospin, and S (I) is the total spin (isospin) 
of N or A. Because the intermediate baryon states are restricted to the lowest mode, the 
self-energy contribution, E^t, has the same form as the HF result except for the spin-isospin 
matrix element: [o^] s / = 27/4 for both N and A. 

The pion-induced baryon self-energies should reproduce the correct, leading non-analytic 
(LNA) behavior of chiral perturbation theory (%PT) 12j. The LNA contribution is as- 
sociated with the infrared behavior of the baryon self-energy, and hence, for example, 
Eq.flED gives -t>g\m\ / '(32tt f 2 ) x (30/25,6/25) for the (N, A) in the infrared limit. In 
contrast, E^^ gives —3g\ml/(327i-f 2 ) x 27/25 for both N and A. Thus, the total amount 
is —3g\ml/(32nf 2 ) x (57/25,33/25) for the (N, A), which is precisely the leading-order 
correction given by large N c x?T. However, because of the closure approximation taken for 
the intermediate states, the term of m^.ln(m„-) does not appear in the present calculation. 

The contribution from the ps (gluon) interaction, E^ Fqs j-(Eg HFqs j-), is also calculated in 
the similar manner. 3 Note that there exists a nonvanishing, interference (sv) contribution 



The energy shift due to the Coulomb propagator vanishes because of the color charge neutrality [13 1 



6 



between the ps and pv couplings, E s qHFqs j. The pv, ps and gluon corrections lower the baryon 
mass, while the interference contribution increases it but its magnitude is Oimj f 2 R) and 
thus small. 

Each correction is a function of the bag radius R, and, for example, E v qliF diverges like 
~ —1/R? as R — > 0. Thus, the bag collapses as R — > 0. Because the pion has a finite size, 
the effect of the qq substructure is essentially important when the bag radius is very small. 



In Ref. 



14j, a phenomenological, non-local interaction was studied to settle this collapse at 



R ~ 0. The effect of the qq substructure of pion can eventually be described by a form 
factor at the vertex of the quark-pion interaction. When the charge radius of the pion is 



about 0.56 fm 



15|, the form factor is estimated as 



Fgn(R) = w ( 15 ) 

q (1 + 1.3 x (b/R) 2 f 2 

with b = 0.46 fm (see Fig. 2 in Ref. 14j). Using this form factor, one can get the finite, 
pion-loop contributions. 

The coupling between a quark and gluon is scale-dependent and the lowest-order coupling 
at momentum transfer Q 2 is a s (Q 2 ) = g 2 /An = 127r/[(33— 2Nf) \n(Q 2 / Aq CD )] with Nf quark 
flavors and Aqcd — 200 MeV. In practice, this behavior can be parametrized in a convenient 
form [16] 

a s (Q 2 ) = a k e~ Q2 ^ = a, e" 4Q2 + 0.25 e' Q " + 0.15 e^ 10 + 0.2 e ~ Q2 ' im \ (16) 
k 

where Q 2 in GeV 2 and the parameters, a k and 7&, except a\ are constrained to follow the 
behavior of a s (Q 2 ). We treat a± as a parameter. The form (|T6|) is convenient, because it is 
easily transformed into the form in coordinate space 



aJR) = Y=£ e~ x2 dx. (17) 



Here we assume that this gives the coupling constant at the scale of the bag radius, R, and 
that the energy shift due to the gluon exchange is given by replacing g 2 with Aira s (R) in 
E 9 qHFqsf (R). Note that a s (R) -> as R -> 0. 

Now we are in a position to present the numerical result 
space. The mass is given by a sum of the usual bag energy 



'or the N or A mass in free 
and the corrections due 
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to the pion and gluon exchanges. We fix the current quark mass m = 5 MeV, because 
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FIG. 1: Energy shift due to the HF contribution (Rn = 0.8 fm). The dot-dashed curve presents 
the gluon contribution, which approaches a constant as R — > 0. The solid (dashed) [dotted] curve 
is for the ps (pv) [sv] contribution from the pion exchange. 

the dependence of the baryon mass on m is very weak. There are four parameters: B, zn, 
za and a\ (in a s ). Since we can expect that the usual z parameter for the N is not much 
different from that for the A, we choose Zq = zn = z&. Then, the bag constant, B, and z 
are determined so as to fit the free nucleon mass, Mn(— 939 MeV), with its radius Rn = 0.6 
or 0.8 fm. The remaining parameter, a 1; is fixed so as to yield the correct mass difference 
between Mn and Ma(= 1232 MeV) together with the pion-cloud contribution. We then 
find B 1 ' 4 = 231.8 (183.7) MeV, z = 2.46 (1.17) and a x = 5.01 (6.08) for R N = 0.6 (0.8) fm. 



TABLE I: Energy shift (in MeV) due to the pion or gluon exchange. 





R(fm) 


T-iPV+pS+SV 

h qHF 




E 9 

qHF 


E 9 . 
qsf 


N 


0.6 


-89.5 


-80.5 


-111.1 


-333.3 


A 


0.611 


-17.7 


-79.7 


110.7 


-332.0 


N 


0.8 


-68.3 


-61.5 


-120.0 


-360.1 


A 


0.823 


-13.7 


-59.1 


119.4 


-358.1 



In Table HJ we present the energy corrections. The N-A mass difference mainly comes 
from the gluon-exchange HF contribution. Note that the N-A mass difference due to the 
pion cloud is about 60 MeV (for R N = 0.8 fm), and that is near the upper limit allowed 
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> 



c 



20 - 
- 

-20- 
-40 - 
-60 - 
-80- 
-100 - 
-120- 
-140- 
-160- 
-180- 
-200- 



_1 i I i l_ 



k ps HF 

/ • 7i pv HF 

.■ tc sv HF 

; gluon HF 

— | . 1 1 1 1 1 . 1 . 1 1 1 1 1 . 1 r- 

0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2 

i?[fm] 



from lattice QCD constraints [17]]. In Fig. [T], for example, we show the HF energy due to the 
pion or gluon exchange as a function of the bag radius. As expected, the ps contribution is 
quite small. The interference is also small and its sign is positive. Because of Eqs.( ll5p and 
(|17p . the energy shift is finite everywhere and thus the total energy for the N or A mass has 
(global) one minimum at a certain R. 

To describe a nuclear matter, we need the intermediate attractive and short-range re- 
pulsive nuclear forces. As the QMC model is based on the one-boson-exchange (OBE) 
picture pj, it is achieved by introducing the a and uj mesons. 4 However, the present a 
meson is chirally singlet and not the chiral partner of the ir meson. This a represents, in 
some way, the exchange of two pions in the iso-scalar channel [rsl ]. 

Now let us start from the following Lagrangian density for the "chiral quark-meson cou- 
pling (CQMC) model": C C qmc = Ccbm + £<tw, where 

C aul = i> \gla - 0« To w] ip 9 V - ^m 2 y + ^m> 2 , (18) 

with g^(g^) the o"(u;)-quark coupling constant, m cr (m aJ ) the meson mass and c{ui) the mean- 
field value of the a (u) meson. 

In an iso-symmetric nuclear matter, the total energy per nucleon is given by 

E tot = ^ 3 / " dkJk? + Mff + 3g*u + \{mla 2 - mlu 2 ), (19) 

with ps = 2kp/3n 2 and the effective nucleon mass. The attractive force due to the a 
changes the quark mass in matter as m* = m — g q a o (in* the effective quark mass), which 
modifies the quark wave function. This modification generates the effective nucleon mass, 
M^r, in matter. Because the change of the quark wave function varies the source of the a field, 
we have to solve the coupled, nonlinear equations for the nuclear matter self-consistently 
(for details, see Refs.jl]). 

The numerical result for the nuclear matter (with m a = 550 MeV, m u = 783 MeV 
and Rn = 0.6 or 0.8 fm) is presented in Table [III The cr-N and o>-N coupling constants, 
g a {= 3g%S(a = 0)) and g u {= 3g q ), are determined to fit the nuclear saturation condition 
(—15.7 MeV) at normal nuclear density po(= 0-15 fm~ 3 ). Here S(o~) is the scalar density 
calculated by the quark wave function [l(|. 



4 It should be noticed that the OBE model is still the most economical and quantitative phenomenology 
for describing the nuclear force [2|. 
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TABLE II: Coupling constants and calculated properties for symmetric nuclear matter at pq. The 
last three columns show the relative changes (from their values at zero density) of the bag radius, 
the lowest eigenvalue and the root-mean-square (rms) radius of the nucleon calculated with the 
quark wave function. The nucleon mass and the nuclear incompressibility, K, are in MeV. 



R N (fm) 




fl£/4* 


M* 
N 


K 


SR* N /R N 


5x*/x 


5r*/r 


0.6 


6.11 


10.20 


632 


362 


0.01 


-0.18 


0.05 


0.8 


4.93 


9.59 


647 


365 


0.02 


-0.22 


0.07 




FIG. 2: Scalar polarizability in the CQMC model (for i?jy = 0.8 fm). The dashed curve presents 
the scalar polarizability in the usual QMC model, while the solid one is for the CQMC model. The 
dotted (dot-dashed) curve shows the pion (gluon) contribution to the scalar polarizability. 

In the CQMC model, the bag radius is swelled by a few percent at po- The quark 
eigenvalue, x, decreases by about 20%, which leads to the smaller in-medium nucleon mass 
than in the QMC model. Although the rms radius of a nucleon swells by about 6% at po, it 



may still be within the experimental constraint [191 ] . 

In Fig. |2l we present the scalar polarizability in the CQMC model, which is given by 
a sum of the quark scalar density, S(a), and the contributions from the pion and gluon 
exchanges. We find that even in the CQMC model the scalar polarizability decreases with 
increasing p#. Because of this reduction of the scalar polarizability in matter, the present 
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FIG. 3: Variation of qa in matter (for Rn = 0.8 fm). 



model can achieve the nuclear saturation property with a much smaller value of the nuclear 
incompressibility, K, than in QHD. 

In Fig. [3X as an example, we show the axial vector coupling constant, gA, in matter. If 
the pion is massless, the pion pole term in the axial vector current gives just 1/3 of the 



quark core contribution i.e., the MIT-bag-model value [10|, [20j. However, when the pion is 
massive, the pion current does not contribute to gA- Here the quark core contribution is 
thus calculated with the corrections due to the Peierls-Thouless projection method and the 
Lorentz contraction of the bag [20]. We obtain gA = 1-15 in free space. The CQMC predicts 
a reduction of about 14% for at po- Note that the pion field surrounding a nucleon 
becomes also weaker with increasing ps since the pion itself is mainly generated by the pv 
interaction. 

In the present model, of importance is the scalar polarizability. It describes the response 
of a quark to the scalar field in matter and leads to the reduction of the <r-N coupling 
constant l[ [21]. Such response arises from the change of the quark wave function in a 
medium, and it is thus the many-body effect Then, the CQMC predicts that the 
iso-scalar, central force is weakened depending on p%. On the other hand, in EFT, the 
(supposed) a exchange in the OBE model can be well understood by the correlated, two-pion 



contribution 



22 



231 ] . However, recent calculations for matter show that, the N 3 LO potential 



based on EFT produces very deep overbinding at large ps [241 ] . Thus, this reduction of the 
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central force may be favorable. To draw more definite conclusions, however, further studies 
are necessary. 

In summary, we have developed for the first time a chiral version of the quark-meson 
coupling model based on the cloudy bag model, in which the effects of pion cloud and gluon 
exchange are included self-consistently. The model can describe a symmetric nuclear matter 
reasonably well. We have also shown that decreases with increasing ps- This implies 
that the iso-scalar, central nuclear force should be weakened in matter. At 0(1/ f^), the 
CBM Lagrangian automatically provides the Weinberg- Tomozawa (WT) term, which is a 



new source of gA [20 



251 ]. Thus, in the future, it is desirable to perform a self-consistent 



calculation up to 0(1/ f^) including the WT term. It is also very interesting to compare the 
CQMC with the EFT approach to investigate the internal structure change of the nucleon 
in medium. 
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